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Abstract. The Kondo lattice model has been analyzed in the presence of a random inter-site interaction 
among localized spins with non zero mean Jo and standard deviation J. Following the same framework 
previously introduced by us, the problem is formulated in the path integral formalism where the spin 
operators are expressed as bilinear combinations of Grassmann fields. The static approximation and the 
replica symmetry ansatz have allowed us to solve the problem at a mean field level. The resulting phase 
diagram displays several phase transitions among a ferromagnetically ordered region, a spin glass one, a 
mixed phase and a Kondo state depending on Jo, J and its relation with the Kondo interaction coupling 
Jk- These results could be used to address part of the experimental data for the CeNi\- x Cu x compound, 
when x < 0.8. 

PACS. 64.60.Cn Order-disorder transformations; statistical mechanics of model systems - 75.10. Nr Spin- 
glass and other random models - 75. 30. Mb Valence fluctuation, Kondo lattice, and heavy-fermion phe- 
nomena 
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1 Introduction 

The magnetism in strongly correlated f-electron systems 
has become a source of great interest due to the physics 
involved Q like, for instance, quantum phase transitions 

[ and Non- Fermi liquid behavior The anti- ferromagnetic 
s-f exchange coupling of conduction electrons to local- 
ized spins can be responsible for the competition between 

, the Kondo effect, that reduces the localized magnetic mo- 
ments, and the RKKY interaction among magnetic impu- 
rities which, in turn, may give rise to magnetic long range 

■ order. 

Recently, an experimental magnetic phase diagram of 
the Kondo CeNi\- x Cu x compound has been proposed 
H showing the existence of a spin glass like state. In the 
CeCu limit, the negative magnetic interaction is domi- 
nant enough to produce an anti-ferromagnetic long range 
order with no indications of the Kondo effect. When Cu 
is substituted by Ni, there is a phase transition around 
X — 0.8 from the antiferromagnetic (AF) to a ferromag- 
netic (FM) ordering, which Anally disappears at roughly 
x = 0.2; the Curie temperature is roughly equal to IK 
and is slowly decreasing down to x — 0.4 and then dis- 
appears at x = 0.2. Above the ferromagnetic phase, a 
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spin-glass (SG) phase was identified by magnetic suscep- 
tibility measurements and the SG transition temperature 
increases from 2 to 6K for x varying from 0.7 to 0.2. For 
x < 0.2 a Kondo behaviour has been proposed, and fi- 
nally CeNi is an intermediate valence compound. Thus, 
at very low temperatures, the phase sequence FM-SG- 
Kondo has been observed with decreasing x and in the 
range 0.7 — 0.2 for x, the sequence FM-SG is obtained 
with increasing temperature. It is quite rare to observe a 
ferromagnetic phase in Cerium Kondo compounds, while 
antiferromagnetic phases are often observed and for exam- 
ple the sequence of SG-AF- Kondo transitions is obtained 
with increasing x in Ce2Aui- x Co x Si3 alloys 

Quite recently, a model has been introduced || to 
study the interplay between spin glass ordering and a 
Kondo state. This model is based on the previously in- 
troduced Kondo lattice model |(| with an intrasite s-f ex- 
change interaction and an intersite long range random in- 
teraction of zero mean that couples the localized spins. 
The use of the static approximation and the replica sym- 
metry ansatz has made possible to solve the problem at a 
mean field level. This fermionic problem is formulated by 
representing the spin operators as a bilinear combination 
of Grassmann fields and the partition function is found 
through the functional integral formalism p,p|,p|,p!0||. The 
results are shown in a phase diagram of T / J versus Jk/ J 
where T is the temperature, Jk is the intrasite Kondo 
exchange interaction and J is the standard deviation of 
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the random inter-site interaction. For high temperatures 
and small values of Jk, a paramagnetic phase is found. 
In this situation, if the temperature is decreased a second 
order phase transition to a spin glass phase appears at X/. 
The model shows a transition line J K (T) separating the 
paramagnetic and the spin glass phases from the Kondo 
phase. 

In the present work, the model mentioned in the pre- 
vious paragraph has been extended in order to include 
the proper elements that produce also a ferromagnetic or- 
dering by taking the mean random interaction Jo to be 
different from zero. Therefore, the magnetization can be 
introduced in addition to the other order parameters and 
solved coupled to them. 

From this procedure a quite non-trivial phase diagram 
is obtained which contains ferromagnetism, a mixed phase 
]l2|,[l3| (ferromagnetism and spin glass), a spin glass phase 
and a Kondo state. For instance, one of the achievements 
of the present work is the finding of a mixed phase whose 
existence should not be discarded in the magnetic mea- 
surement of CeNii^ x Cu x , as mentioned in Ref. || 

This paper is organized as follows. In section II we 
present the model and its development in order to obtain 
the free energy and the saddle point coupled equations for 
the order parameters. The phase diagram of the temper- 
ature T/J versus Jk/J is shown for several values of Jo. 
The Almeida- Thouless line is also calculated. Discussions 
and concluding remarks are presented in the last section. 



2 The model and results 

The model considered in this work was introduced before 
in Ref. || to study spin glass ordering in a Kondo lattice 
compound so the Hamiltonian is 



The /if (pL c ) are the chemical potential for the localized 
(conduction) band. The energy eo is referred to /// while 
e& is referred to fi c . 

The partition function is expressed in terms of func- 
tional integrals using anticommuting Grassmann variables 
(fiia- (r) and ipia (t) associated with the conduction and the 
localized electrons respectively. Therefore, 



Z= j D(^)I%V)(sq> M dr[L a (r^) + 

Lo&*,<p)+L SG + L K }\ (3) 



where 
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in the static approximation |t|,^,[|[io), it is possible to 
solve the problem in a mean field theory where the Kondo 
state is described by the complex order parameters [pl|6|: 
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where Jk > and the sum runs over N lattice sites. In 
the present case the random intersite interaction Jjj in 
the Hamiltonian is infinite ranged with a Gaussian dis- 
tribution where < Jy- >= 2J /N and < J?j >= 8J 2 /N. 
This particular scaling compensates the factors 1/2 that 
originate in the definition of the operators S z in equation 
@ and also in changing from sum over sites to sum over 
bonds. 

The spin variables S^' (s^ S% are bilinear com- 
binations of the creation and destruction operators || for 
localized (conduction) fermions j\ al fi a (dj^dj^) with the 
spin projection a =| or [: 



4i = 4t d n 

S ci = 4jAT 



\[fy^-fihi\ (2) 



i.tu 
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Following the treatment for the Kondo part in the par- 
tition function as introduced in Ref. where it was as- 
sumed that A* ~ A* (A CT sj A), we show in the Appendix 
that first the conduction electron degrees of freedom may 
be integrated out to give 
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where 



Z eff = J exp EE^H^M^H 
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SG 
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Z® is the partition function of the free conduction elec- 
trons, 



(8) 



while 7 J ; J 1 = iuSij — PUj is the inverse d-electron Green's 
function and (3 = 1 /T is the inverse temperature. 
The free energy is given by the replica method 

pF = 2(3J K X*X - Urn J-(((^ // (J, J .))) co - 1) (9) 

and the averaged replicated partition function can be lin- 
earized by means of the usual Hubbard-Stratonovich 
transformation. Therefore, 



n a f3dq al 3 J II a dm a exp 



-N x 



P 2 J 2 



E2 . PJo sr^ 
i«0 + ^rz^ 



a0 




A(q a p,m a ) (10) 



with a — 0, 1, ..n being the replica index and 
A(q al3 ,m a ) = J D{ip* a ip a ) expi ^ 

C« HVW M + P Jo £ 2 ^>" 



/? 2 J 2 £ 4Sf Sf <^ > • (11) 
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A more detailed derivation of equations (10) and ( |ll| ) is 
given in the Appendix. 

This problem is analysed within the replica symmetric 
ansatz where q a ^/3 = q is the spin glass order parameter, 
m a = m is the magnetization and q aa — q + x, (x = 4) 
with x being the static susceptibilty. The sum over replica 
indices also gives quadratic terms which can be linearized 
again by introducing new auxiliary fields in equation ( |l0| ) : 

A(q a p,m a ) = J Dzj J D(^*^p) x 

act y-^n ) ~r 

I i,ja 7 uj a 



( exp {^f + PJ^H* 2S« 

J [ <*j 



(12) 



The functional integral in equation (|12|) can be per- 
formed and the saddle point solution for the free energy 
is given by 



1 



(3F = 2(3 J K \ 2 + -P 2 J 2 {x 2 + 2qx) 



PJ 



lim — — 

n->o Nn 



n,,Dzi x 



J II a>i Dtf exp 



J> [det Gr»] 



- 1 



(13) 



where in the previous equation we introduced the inverse 
Green's function 



G7 j l(cj)=gT j 1 (u J )-6 ij ah(z i ,g) 
with an effective field 



h(z h £?) = pJ m + pJy/2qzi + 0Jy/2x& 



(14) 



(15) 



A problem is presented by the calculation of the 
Green's function Gij a (uj) in equation (jlj) where there is a 
random Gaussian field h(zi, £f ) = hi a applied at every site 
i of n replicated lattices with N sites. The decoupling used 
at this point is the same as for Ref. i.e., the original 
Green's function Gij a {u) is replaced by a Green's function 
r^ua{^) where there is a uniform field h ia applied in ev- 
ery site fj,, v of a ficticious Kondo lattice. Therefore, going 
to the reciprocal space and assuming a constant density 
of states for the conduction electron band p(e) = for 
—D<e<D, the sum over Matsubara's frequencies oj 
can be performed in equation (|l^) and the resulting free 
energy is 

0F = 2PJ K \ 2 + U 2 J 2 {x 2 + 2xq) + ^m 2 - 



Dz In 
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D£ e 
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where 



+0D 



P D J-0D 



S(£,x) = cosh(/i+) + cosh(VZi), 

A = (h-) 2 + (pj k X) 2 and ht = 

The saddle point equations for the ord er p arameters q, m, 
X and A can be found from equation (|l6|). 

The limit of stability for the order parameters solutions 
with replica symmetry is achieved if the Almeida-Thouless 
eigenvalue becomes negative: 



\ AT = \-2P 2 J 2 \ Dzl J -°°* — 



where Dx = —7=e x2 / 2 dx and 

V27T 



u) n a.(7— ± 



(17) 
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Fig. 1. Cut in the phase transition space transversal to the 
Jk/J axis for Jk= 2, J = 0.5 and D — 10. The Kondo state 
is not turned on yet and the transitions among the ferromag- 
netism and the spin glass phases depend on the value of J . The 
dashed line shows the transition from the paramagnetic phase 
to the ferromagnetic and the spin glass phases. The dotted 
line is the Almeida- Thouless (AT) line which, for lower values 
of Jo, coincides with the paramagnetic - spin glass transition 
line (horizontal dashed line). The dot-dashed and the AT line 
delimit the mixed phase between the ferromagnetic and spin 
glass phases. 
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3 Discussion 

We have studied in this work a Kondo lattice model where 
the localized moments interact through a random inter- 
site interaction which has an average different from zero. 
The static approximation and replica symmetry ansatz 
lead to a mean field solution for the problem. The result- 
ing coupled saddle point equations for the order parame- 
ters produce solutions which give a Kondo state and mag- 
netic ordering like ferromagnetism, spin glass and a mixed 



phase. In principle, we would be able to build up transition 
surfaces among those phases in a space T / J (temperature) 
versus Jo/J (the inter-site interaction average) and Jk/J 
(the Kondo coupling) where J (the inter-site interaction 
standard deviation) is kept constant. These parameters 
Jk, Jo, J and the temperature are the set of energy scales 
in the present model. The conduction electrons bandwidth 
D is kept constant. 

The result shown in Figure [l] represents a cut in the 
cited space transversal to the Jk/J axis. For this situ- 
ation and using Jk/J= 4, the Kondo state is still not 
turned on (it means A = 0). The obtained phase diagram 
for this fermionic model resembles the classical one [ pT| , p^ | 
(it depends basically on Jo and its relation to J) except 
that the numerical values of the transition temperatures 
are smaller. If Jo < 1.46 J, for decreasing temperature, 
there is a second order transition from a paramagnetic to 
a spin glass phase (m = 0,? ^ 0). For that region, the 
Almeida-Thouless (AT) line coincides with the transition 
line. Howewer, for Jo > 1.46 J, for decreasing tempera- 
ture, the model shows a transition from a paramagnetic 
to a ferromagnetic phase (m / 0, ^ / 0). The AT line is 
located at higher temperature than the calculated replica 
symmetry line transition T*(Jo) between the ferromag- 
netic and the spin glass phases. Therefore, this fermionic 
model shows a transition from a ferromagnetic to a replica 
symmetry breaking spin glass phase with a large number 
of degenerate states but still with non-zero spontaneous 
magnetization which is called a mixed phase p2 ,13. 

Figure |^ shows the cut in the phase space transversal 
to the Jo/J axis. For values of Jo/J close to zero (see Fig- 
ure ||a) , the phase diagram resembles the scenario already 
found in Ref. [||, that is a paramagnetic phase at high 
temperatures, a spin glass phase below the freezing tem- 
perature Tf and a line Jk =J k (T) separating both phases 
from a Kondo state. For that situation, the AT line is at 
T f and follows the line J K =J K (T). 

As the value of Jo/J is increased (for small Jk/J) 
(see Figures ||b, ^c and ||d), the phase diagram starts to 
show the presence of a ferromagnetic phase which has a 
transition temperature T c (Jo) increasing with Jo, the AT 
line and the calculated replica symmetric line T*(Jo) de- 
creasing with Jo. Nevertheless the AT line is always above 
T*(Jo). In that scenario, for decreasing temperature, first 
a transition from paramagnetic to a ferromagnetic phase 
appears followed by a transition from the ferromagnetic to 
a mixed phase . This behavior is reminiscent of that one 
described for the cut transversal to Jk/J (Figure 0). For 
some value of Jo / J, the mixed phase finally disappears 
and that region of the phase diagram is totally occupied 
by the ferromagnetic phase. For larger values of Jk/J, the 
phase diagram goes to a Kondo state where the transition 
line J£-(T) does not depend on J . 

A remark should be made about the transition line be- 
tween the spin glass phase and the Kondo state. At low 
temperatures, this is a first order transition line and so 
multiple possible solutions for the order parameters can 
be found. Nonetheless, the actual stable solutions can be 
obtained from the minimization of the free energy. In the 
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Fig. 2. Cut in the phase transition space transversal to the Jo/ J axis for several values of Jo, J = 0.5 and D = 10. The solid 
line shows the thermodynamically stable transition from the Kondo phase to the other ones. For lower temperatures, a hatched 
region in panel (a) delineates a multiple solution region where the solid line to the right corresponds to the thermodynamically 
stable solution. The horizontal dashed, dotted and dot-dashed lines have the same meaning as in Figure |l[ The AT line follows 
the horizontal dotted line up to the Kondo transition point. Beyond that, for larger values of Jk, the AT line follows the 
transition line from the Kondo phase and the mixed and spin glass phases (solid line, lower temperatures). One can notice 
that, as Jo increases, a ferromagnetic and a mixed phase start to appear and, for some value of Jo, the spin glass phase finally 
disappears. 



case of Jo/J< 1.46 the hatched region in Figure ||a dis- 
plays where these multiple solutions occur. By comput- 
ing the free energy we have found the thermodynamically 
stable solutions. The solid line to the right of the hatched 
region in Figure ||a corresponds to these solutions and the 
hatched region itself corresponds to metastable solutions. 
Such a carefull analysis can be considered an improvement 
with respect to our previous work || where such a discus- 
sion had not been done. Nevertheless the previous SG- 
Kondo state transition line shown there is approximately 
the same as the one presented here. Thus, the hatched 
regions corresponding to the one displayed in Figure 
have not been presented in Figures ^b-d. 

4 Conclusions 

In this work it has been studied a Kondo lattice model 
in the presence of a random inter-site interaction which 
produces paramagnetism, ferromagnetism, a spin glass 
phase, a mixed phase and a region where the magnetic 
moments of the localized electrons are supressed by the 



screening of the conduction ones (Kondo state). The 
model has four energy scales: the temperature, Jo (the 
average inter-site interaction), J (the inter-site interac- 
tion variance) and Jk (the Kondo coupling). As a re- 
sult one has a three dimensional phase diagram with axes 
Jo/ J, Jk/J and temperature. Some cuts of this diagram 
transversal to the Jk / J and Jo/ J, planes are shown in the 
Figures [j] and ^. The position of the transition line sep- 
arating the Kondo state from the magnetic phases is not 
affected by Jo- This energy scale is basically responsible 
for locating several magnetic orderings along the temper- 
ature range. 

One can try to address the experimental phase dia- 
gram found in Ref. || for the alloys CeNii_ x Cu Xl but 
theoretically if we vary only Jk with x, we have found 
a ferromagnetic phase above the spin glass, in disagree- 
ment with the experimental result. However, the equiva- 
lence between their experimental phase diagram and ours 
(see Figures |l| and |^) is not so straightforward since the 
Ni content would have to be associated to both Jo and 
Jk ■ This could be an indication that the ergodicity break- 
ing mechanism for the formation of magnetic phases like 
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spin glass and ferromagnetism in CeNi\- x Cu x is far more 
complicated than the modelling by a random inter-site in- 
teraction can address. Although recent investigations on 
the ferromagnetic transverse Ising spin glass suggest also 
the existence of a spin glass transition below the Curie 
temperature [H , it is plausible that this be a characteris- 
tic of the Sherrington-Kirkpatrick model with a high de- 
gree of frustration. Less frustrated spin glass models M] 
may sustain spin glass order above the Curie tempera- 
ture and they can be more indicated for the study of the 
CeNi\^ x Cu x compounds. 

To conclude, in Ref. Q a Kondo lattice model with 
strong frustation (at mean field level) has been solved 
showing the existence of a SG and a Kondo state de- 
pending on Jk/J (as defined in Section ||). These re- 
sults could address part of the experimental phase di- 
agram of CeNii^ x Cu x j|. The purpose of the present 
work has been to examine a wider and more complex re- 
gion of this experimental phase diagram which includes 
ferromagnetism. Therefore, we have improved our previ- 
ous work by choosing a non-zero average of the random 
coupling Jo. From this approach we have been able to 
generate a quite non-trivial phase diagram with a spin 
glass phase, ferromagnetism, a Kondo state and a mixed 
phase (spin glass and ferromagnetism). Nevertheless the 
calculated spin glass freezing temperature is lower than 
the Curie temperature in contrast with some experimen- 
tal findings ||. However, as pointed out in Ref. a 
mixed phase can not be discarded as a possible expla- 
nation for the magnetic measurements. The calculations 
with the ferromagnetic phase has also shown an improve- 
ment with respect to our previous work Q regarding the 
actual location of the SG-Kondo first order transition line. 
The present approach might also explain the frustration in 
antiferromagnetic Kondo systems like CeiAu\- x Co x Sij, 
alloys. This work is now on progress. 
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Appendix 

We outline here briefly the method used in Rcf.j5[. In or- 
der to obtain equations (]l(]) and ( |Tl| ) of Section p we must 
first use the static approximation in the Fourier transform 
of Lk in equation (|4|) and introduce the Kondo order pa- 
rameter by means of the identity 



cxp ^ J Lk<It 



IlcrdX^dXcync X 



(w)ip ja (w)) X 



cxp 



(w-Eiw^w) x 

{WV^ + aJAt]} (18) 



and using the integral representation for the 5 function we 
obtain, after some algebra, 

Z = J n a d\ld\ a exp |-iV/3J K A - A -| ( 19 ) 



where 



Zstat = \ D^)D{ip*ip) X 



exp I J dr [L (ip*, VO + Lo(tp*, <p) + L SG ] \ x 



exp < 13 J K ^ 



(20) 



The mean field approximation adopted here is based 
on two assumptions: first, fluctuations in time are ignored 
(static approximation); second, fluctuations in space are 
also ignored in the definition of the order parameters. Both 
assumptions lead us to a quadratic form in the Grassmann 
variables tp and ip* in equation ( pp| ) and shifting the repre- 
sentation to Matsubara's frequencies we can perform the 
functional integrals to obtain equation (j(^) , where now the 
order parameters AJ., X a are taken at their saddle point 
value. We also have 



drLsa 



7 / Ji J S h S h 



(21) 



where 



s h = \ E [$r H^t H - i>ti M<M H] ■ (22) 



Hence, in order to get the configurational average over 
the random coupling J^-, we use a gaussian distribution 
with average and variance given in Section ^. So 



where 



D(^V)exp{Af +A r s f ic ] (23) 



A = EE^^^ii 1 ^^^^"' 



S. G. Magalhaes et al.: Spin glass and Ferromagnetism in Kondo lattice compounds 



7 



and 



Areplic 1 J /? J \ 

a sg \ in 2^ 



a/3 



-i 2 



In the previous equation, g^ 1 (w„) has been defined in Sec- 
tion ^| The static approximation is also used to write the 
S" in terms of Grassmann fields as in equation ( p2| ) and 
the resulting equation can be linearized by standard pro- 
cedures 0, introducing the order parameters q a p and 
m Q which gives equations (|ic| ) and (pi). 
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